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Abstract. A generalization of the Gell-Mann-Low theorem is applied to the antimuon-electron 
system. The bound state spectrum is extracted numerically. As a result, fine and hyperfine structure 
are reproduced correctly near the nonrelativistic limit (and for arbitrary masses). We compare the 
spectrum for the relativistic value a = 0.3 with corresponding calculations in light-front quantiza- 
tion. 
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More than thirty years have passed since quantum chromodynamics (QCD) has been 
formulated. Overwhelming evidence has accumulated during this time that QCD gives 
a complete description of hadronic physics at every energy scale currently accessible 
to experiment. However, the status of the theory is not at all satisfactory: we are still 
searching for a deeper understanding of the interaction between quarks and gluons 
that arises from the QCD Lagrangian beyond the high momentum transfer limit where 
perturbation theory is applicable. Parts of the physical spectrum that is generated by the 
theory can be determined through lengthy numerical calculations on space-time lattices 
which, however, do not provide major insights into the underlying dynamics. Ideally, 
one would like to have an analytical description of the interaction between quarks and 
gluons derived from the QCD Lagrangian, typically in the form of an approximation 
and a systematic procedure to incorporate corrections to the latter. No such description 
is available to date. Moreover, once it is obtained, we still face the problem of calculating 
the bound states resulting from this interaction, i.e., the physical hadrons. 

The present contribution deals with this second step in connecting the theory in a 
transparent way with the phenomenology. Without an appropriate description of the fun- 
damental interaction at hand, we consider a simpler theory, namely, quantum electrody- 
namics (QED), which still has some features in common with QCD. In order to simulate 
a relativistic situation like the one that prevails in the lighter hadrons, we artificially con- 
sider larger coupling constants. E.g., a = 0.3 was used before in light-front calculations 
we will compare with later on. For concreteness, we will consider muonium, a bound 
state of an antimuon with an electron, although we will allow for an arbitrary antimuon 
mass. In the special case of equal masses of electron and "antimuon", the system is 
similar to positronium if we disregard all (virtual and real) annihilation processes there. 

Our approach to relativistic bound states is an application of a generalization of the 
Gell-Mann-Low theorem [1] to the subspace of Fock space that contains all states of 
one electron and one antimuon (and no photon). As in the earlier applications of the 



same formalism to the Wick-Cutkosky model and Yukawa theory [2, 3], the effective 
Hamiltonian generated by the generalized Gell-Mann-Low theorem contains the rela- 
tivistic kinetic energies of the constituents and an effective potential. To lowest order in 
a perturbative expansion, the matrix elements of the effective potential read, in Coulomb 
gauge. 
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Here, |Pa,'";Ps,>s') symbolizes the state of an electron with 3-momentum pa and spin 
orientation r (in a spinor basis yet to be specified) and an antimuon with 3-momentum 
Pb and spin orientation s. We use the shorthands E^^ = (m^ + p^)^/^ and E^^ = 
(m| + p|) for the kinetic energies. For convenience, we have introduced the charge- 
conjugate Dirac spinors ub{'Pb,s) for the antimuon, while ua{'PA:S) represents the elec- 
tron spinors. The spatially transverse photon polarization vectors ej-^\ii) satisfy the re- 



lation Lj^i £^^) (k) e)"^'"(k) = 55(k) = 5ij - k,kj (where k = k/|k|). 

As for the interpretation of the effective potential (1), the second line stems from the 
instantaneous Coulomb potential, easily identified by the momentum dependence in the 
denominator (the Fourier transform of the spatial Coulomb potential), and multiplied 
with the charge densities of the Dirac currents. The following lines are the result of 
transverse photon exchange, the more complicated denominator indicating a retarded 
interaction, and the Dirac currents being contracted with the corresponding photon 
polarization vectors. 

The delta function in Eq. (1) shows that total 3-momentum is conserved by the 
effective interaction, and in the following we will consider the center-of-mass system 
(c.m.s.) Pa + Pb = Pa + Ps = 0- order to simplify the diagonalization of the effective 
Hamiltonian, we express the Dirac spinors in terms of Pauli spinors (using the Dirac- 
Pauli representation) to find for the effective Schrodinger equation in the c.m.s.. 
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Here, p = p^ = — ps, the spinorial wave function ^(p) is defined as 0(pA)(2;r)-'5(pA + 
Pe) = Lr,i(PA,^;Pfi,>^|0) [Xr®Xs], and cr^ (Cs) is understood to act on the Pauli spinor 
Xr iXs) only. E' is the difference between the energy of the bound state and the vac- 
uum energy. Of the full state |0) in Fock space (with zero total momentum), only its 
projection to one-electron-one-antimuon states |pA,r;pB,s) appears. The effect of its 
components in other Fock space sectors is taken care of implicitly by the effective po- 
tential. 

In order to solve the Schrodinger equation (2), we take into account its rotational 
symmetry. Eigenstates of total angular momentum J can be constructed as usual by 
adding relative orbital angular momentum L and total spin S. For convenience, we will 
label the eigenstates by the "relative parity" Tz' defined through (— 1)^ = :;r'(— 1)''. Since 
S = 0, 1 , for given / the sector k' — contains the states with L — J and 5 = or 5 = 1 , 
while for %' = —1 we can have L = J— \ or L = J +\, with 5 = 1 in both cases. In any 
sector , the two different possible (LS) -states will mix, except in the following cases: 

(i) for 7 = 0, r' = 0+ is only realized by (L = 0, S = 0) , and 0" only by (L = 1 , 5 = 1 ) ; 

(ii) in the case of equal masses, the Hamiltonian acquires an additional symmetry under 
the exchange of particles A and B; as a result, S becomes a good quantum number and 
there is no mixing in the {Tt' = -|-l)-sector; (iii) in the one-body limit where one of the 
masses goes to infinity, the spin of the heavy particle decouples from the dynamics; as 
a result, every two states are degenerate in this limit and L becomes a good quantum 
number [no mixing in the (tt' = — l)-sector]. 

After explicitly carrying out the contractions of the spatial indices in the transverse 
photon exchange part, the formulae derived before for the application to Yukawa theory 
[3] can be used to determine the result of the application of the terms containing the 
Pauli matrices in Eq. (2) to the total angular momentum eigenstates. On the other 
hand, the application of the factors containing |p — p'| on orbital angular momentum 
eigenstates proceeds through the partial wave decomposition of the former. The partial 
wave decomposition of the (Fourier transformed) Coulomb potential is well-known, the 
one for the 5/y-part of the transverse photon exchange has been calculated in Ref. [3]. 
The partial waves of the kikj-part of the transverse photon exchange are given by 
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This integral diverges like {p — p')^^ for p' p which would lead to a divergence in 
the p'-integral. These divergences, of course, are spurious and cancel in pairs. However, 
for the numerical calculation, we have to extract the divergent parts and perform the 
cancellations analytically. Fortunately, the extraction of the divergencies is simple: they 
occur at cos = 1, and since Pl(1) = 1, we can define "reduced" Legendre polynomials 
Pi {cos e) through 

Pl(cos0)-1 = (cos0-l)if (cos0). (4) 

Separating the one on the l.h.s. of Eq. (4) under the integrals (3), the remainder of the 
integrals is logarithmically divergent for p' ^ p (as are the other partial waves), and the 
following ^^'-integration is convergent. The divergent parts in Eq. (3) originating from 
the one in Eq. (4) can be analytically cancelled in pairs, leaving a finite contribution. 

Putting everything together, the potential term in the Schrodinger equation (2) reduces 
to a one-dimensional integral over p' when applied to the total angular momentum 
eigenstates. For every sector , two such equations are coupled (with the exceptions 
mentioned above). The explicit expressions for the integral kernels (diagonal and off- 
diagonal because of the coupling) are quite lengthy and cannot be reproduced here due 
to lack of space. 

The (coupled) one-dimensional integral equations can be solved by expanding the 

wave function in an appropiate orthonormal basis. After reducing the basis to a finite 
number of elements (40 in our calculations), the integral equations are approximately 
replaced by finite matrix equations. The matrix elements are two-dimensional integrals 
which are calculated numerically (we use a two-dimensional grid of 400 x 800 points). 
Finally, the matrices are numerically diagonalized to give the (approximate) eigenvalues 
and eigenstates of Eq. (2). The results for the lowest energy eigenvalues are plotted 
in Figs. 1 and 2 for equal constituent masses and fine structure constants a < 0.45. 
The binding energies are normalized to jua^, jU being the reduced mass, so that the 
comparison with the nonrelativistic energy eigenvalues jia^/ln^ is immediate. 

We find that for values a < 0. 1, the energy levels are dominated by the nonrelativistic 
values plus the leading relativistic corrections (the leading -order fine and hyperfine 
structure) of order iia^, both in our numerical results and in the perturbative calculations 
of bound-state QED. In this region of small coupling constants, the numerical results 
are in good agreement with the perturbative calculations, apparently only limited by 
the numerical precision. For larger values a > 0.1, higher perturbative orders become 
important and our numerical results deviate in some cases strongly from the lowest-order 
perturbative predictions. 

In Table 1 we compare our results for a = 0.3 with two different calculations in light 
front quantization [4, 5] (we use the data for the Gaussian similarity function in the latter 
paper). In the table, we label the states by the nonrelativistic notation n^^^Lj and also 
indicate the corresponding sectors . There is a clear tendency in our results towards 
more negative energies, i.e., stronger binding, compared to i^(/ia^)-perturbation theory. 
The ordering of the different levels, however, is the same as in perturbation theory. We 
can also see that the difference to perturbation theory in the direction of stronger binding 
is systematically larger for 5-states than for P-states, and also larger for {J = 0)-states 
than for (7 = 1) -states, and smallest for the (7 = 2) -state. For the light-front results, this 
latter tendency is inverted; the So-states have even higher energies than in perturbation 
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FIGURE 1. The binding energy £ = £" — wa — ws as a function of the fine structure constant a = e^/4n 
for the case of equal masses niA = niB- E is normalized to lJ.a^ where [i is the reduced mass (p. = mA/2 
in the present case). Plotted are the lowest energy levels for J'' = 0+ and 1~ corresponding to the 
nonrelativistic principal quantum number n=l. 



TABLE 1. Binding energies £'//ia'^ for equal masses from per- 
turbation theory to ^(;U a"*), from our numerical results, and from 
Refs. [4] and [5] 
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theory. Both light-front calculations are qualitatively similar, only that the binding is 
stronger throughout in the similarity transform approach of Ref. [5]. In conclusion, we 

find very different results with the two different methods for relativistic bound state 
calculations (in the approximations presently considered). We should remark, however. 
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FIGURE 2. As Fig. 1, but for the energy levels corresponding to the nonrelativistic principal quantum 
number n = 2. 




that there is an unphysical logarithmic UV cutoff dependence in the light-front results 
(for the cited values, the cutoff has been set equal to the constituent masses). 
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